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ABSTRACT 

This  report  describes  the  development  of  a  thermal  wave  interferometry  system 
for  the  characterisation  of  film  thickness.  The  system  furnishes  a  measurement 
of  the  relative  phase  between  a  modulated  heat  flux  applied  to  a  thin  film  by  a 
diode  laser  and  the  oscillatory  temperature  response  of  the  film.  The  relative 
phase  is  then  converted  to  a  thickness  estimate  by  means  of  an  inverse  solution 
of  the  heat  diffusion  equation.  A  validation  of  the  facility  was  performed  on 
an  experimental  test  coupon  prepared  with  paint  layers  of  varying  thickness. 
Estimates  of  the  paint  thickness  profile  furnished  by  the  facility  were  compared 
to  and  corroborated  by  independent  measurements  obtained  using  an  eddy 
current  probe  and  a  surface  profilometer.  The  system  is  shown  to  resolve 
variations  in  paint  thickness  of  less  than  1.5  //rri. 
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Development  of  a  Thermal  Wave  Interferometry  System  for 

Thin-Film  Characterisation 

Executive  Summary 

This  report  describes  the  development  of  a  Thermal  Wave  Interferometry  (TWI)  system 
for  the  characterisation  of  film  thickness.  The  capability  addresses  the  need  for  a  method 
to  assess  bondline  uniformity  for  the  compact  multi-parameter  load  evaluation  (CMPLE) 
device,  a  low-cost  surface-mounted  sensor  package  with  an  in-flight  strain  measurement 
capability.  The  device  was  conceived  and  developed  by  the  Defence  Science  and  Technology 
Organisation  (DSTO)  and  successfully  tested  in-flight  on  the  tail  plane  of  an  Australian 
Defence  Force  (ADF)  Caribou  aircraft.  However,  in  validation  testing  of  the  device  it 
was  revealed  that  accurate  strain  readings  require  that  the  device  be  attached  with  an 
adhesive  bondline  that  is  relatively  thin  and  of  uniform  thickness.  In  practice,  this  can  be 
difficult  to  achieve  even  under  nominally  ideal  experimental  conditions,  raising  the  need 
for  a  reliable  method  of  verifying  the  thickness  and  uniformity  of  the  bondline. 

As  an  acutely  sensitive  method  of  measuring  film  thickness,  thermal  wave  interfer¬ 
ometry  has  the  potential  to  provide  an  effective  in  situ  capability  for  adhesive  bondline 
metrology.  As  a  first  step  toward  investigating  that  possibility,  a  laser-based  TWI  facility 
was  developed  and  experimentally  validated.  The  system  furnishes  a  measurement  of  the 
relative  phase  between  a  modulated  heat  flux  applied  by  a  diode  laser  and  the  oscillatory 
temperature  response  of  the  object  under  inspection.  The  phase  measurement  is  converted 
to  a  thickness  estimate  by  means  of  an  inverse  solution  of  the  heat  diffusion  equation. 

A  validation  of  the  facility  was  performed  on  an  experimental  test  coupon  prepared 
with  paint  layers  of  varying  thickness.  Estimates  of  the  paint  thickness  profile  furnished 
by  the  facility  were  compared  to  and  corroborated  by  independent  measurements  obtained 
using  an  eddy  current  probe  and  a  surface  profilometer.  The  system  was  shown  to  resolve 
variations  in  paint  thickness  of  less  than  1.5  /mi.  Although  the  main  purpose  of  the  facility 
is  to  enable  the  development  of  an  in  situ  TWI  capability  for  bondline  metrology,  the 
facility  also  offers  the  potential  to  assist  in  the  evaluation  of  coating  systems  used  in 
military  platforms  for  corrosion  prevention  and  thermal  barrier  applications. 
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1  Introduction 

This  report  describes  the  development  of  a  Thermal  Wave  Interferometry  (TWI)  system 
for  the  characterisation  of  thin  film  properties.  The  need  for  such  a  system  was  uncovered 
during  an  investigation  of  the  compact  mult-parameter  load  evaluation  (CMPLE)  device 
[1],  a  sensor  package  designed  to  provide  a  rapidly  deployable  operational  loads  monitoring 
capability  for  military  aircraft.  The  device  was  recently  applied  in  a  successful  in-flight 
strain  survey  of  the  tail  plane  of  an  ADF  Caribou  aircraft.  To  function  effectively  as  a 
strain  sensor  the  device  needs  to  be  attached  to  a  host  with  an  adhesive  bondline  that 
is  relatively  thin  and  of  uniform  thickness,  otherwise  biased  strain  readings  can  occur. 
However,  in  practice,  a  uniform  bondline  thickness  can  be  difficult  to  achieve  even  under 
nominally  ideal  experimental  conditions,  raising  the  need  for  a  reliable  method  of  verifying 
the  thickness  and  uniformity  of  the  bondline. 

The  installation  procedure  for  the  CMPLE  device  normally  results  in  a  bondline  thick¬ 
ness  of  approximately  50  /jin.  In  general,  acoustic  methods  for  thickness  measurement 
based  on  time  of  flight  are  difficult  to  implement  for  such  short  lengths,  and  acoustic 
resonance  techniques,  although  sensitive,  lack  the  robustness  required  for  this  type  of 
measurement.  Within  that  range  of  thickness,  heat  diffusion  offers  a  more  practical  basis 
for  measurement.  The  use  of  TWI  to  characterise  the  thickness  and  thermal  properties  of 
thin  films  is  well  established  [2,  3,  4],  and  follows  the  pioneering  work  of  Rosencwaig  and 
Gersho  [5]  several  decades  ago.  Thermal  wave  interference  occurs  in  layered  structures  in 
response  to  a  modulated  heat  source  and  describes  a  condition  where  the  properties  of  the 
interface  influence  the  phase  and  amplitude  of  the  temperature  oscillations  at  the  surface. 
In  the  initial  work  of  Rosencwaig  and  Gersho  as  well  as  others  [2] ,  the  surface  temperature 
oscillations  were  resolved  by  means  of  the  photoacoustic  method  which  involves  the  mea¬ 
surement  of  pressure  variations  due  to  thermal  expansion  and  contraction  of  the  boundary 
layer  of  gas  (see  Figure  1).  However,  with  the  widespread  availability  of  low-cost  infrared 
detectors,  photothermal  radiometry  is  now  preferred  over  the  photoacoustic  method  for 
applications  involving  the  measurement  of  thin  film  properties  [3] . 

This  report  describes  the  development  and  validation  of  a  TWI  facility  that  employs 
photothermal  radiometry.  The  main  objective  of  this  facility  is  to  provide  an  experimental 
platform  for  the  development  of  an  in  situ  bondline  thickness  sensor,  however  the  facility 
has  other  important  potential  applications  to  thin-film  material  characterisation,  one  of 
which  is  demonstrated  as  part  of  an  experimental  validation  process  described  in  the 
present  report. 
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backing  (b)  sample  (s)  gas  layer  (g) 


Figure  1:  Layered  structure  used  in  the  work  of  Rosencwaig  and  Gersho. 


2  Theory 

Thermal  wave  interference  is  fundamentally  a  result  of  heat  diffusion  which  is  governed 
by  the  equation, 


--av2r  =  o,  (i) 

where  T  represents  the  change  in  temperature,  t  represents  time  and  a  is  the  material 
thermal  diffusivity.  Rosencwaig  and  Gersho  [5]  show  that  for  the  layered  structure  depicted 
in  Figure  1,  modulated  laser  irradiation  of  the  surface  leads  to  a  surface  temperature 
oscillation  with  the  following  complex  amplitude, 


9 

where, 
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Here,  fd  is  the  optical  absorption  coefficient  for  the  sample,  I o  is  the  incident  light  flux, 
ai  is  the  thermal  diffusion  coefficient  of  the  material  and  l  is  the  sample  thickness.  The 
subscript  i  can  take  the  values  g,  s  or  6  for  the  gas,  sample  and  backing  layers  respectively, 
and  w  denotes  the  modulation  frequency  of  the  incident  light. 

Equation  2  takes  into  account  the  thermal  properties  of  each  layer  within  the  structure 
and  also  the  reflection  and  transmission  of  heat  at  the  layer  boundaries.  Information  about 
the  dimensions  and  thermal  properties  of  a  layered  structure  is  contained  in  the  phase 
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Figure  2:  Semi-infinite  solid  exposed  to  a  surface  heat  flux. 


difference  between  the  source  modulation  and  the  temperature  oscillation  at  the  surface, 
which  is  given  by  the  argument  of  Eqn  2,  i.e., 

</>  =  arg(6>),  (3) 

It  is  important  to  note  that  the  temperature  oscillation  at  the  surface  of  a  semi-infinite 
body  lags  the  heat  source  by  45°.  The  result  underpins  much  of  the  experimental  work 
described  later  in  this  report  so  is  worthwhile  elaborating.  The  result  is  derived  from  a 
harmonic  solution  to  Eqn  1  for  the  semi-infinite  solid  shown  in  Figure  2.  For  consistency 
with  Ronsencwaig  and  Gersho,  the  heat  flux  q(t )  is  given  the  form, 


q(t )  =  /0(1  +  cos(tuf)). 


(4) 


A  harmonic  solution  to  the  Helmholtz  form  of  the  heat  diffusion  equation  can  be 
obtained  by  a  separation  of  variables, 


T  =  9(x)ejui\ 


(5) 


where  9{x)  is  the  amplitude  of  the  temperature  oscillation  as  a  function  of  penetration 
depth.  Substitution  of  Eqn  5  into  Eqn  2  gives, 


P-^0  =  o, 


(6) 


which  has  the  general  solution, 


Q  =  Ae~v^x  +  B. 

Notice  that  6  — »  0  as  x  — »  oo,  so  B  =  0,  reducing  Eqn  7  to, 


(7) 


=  Ae 


(8) 
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The  surface  boundary  condition  involves  a  prescribed  flux,  i.e., 


,  d6 
Q  =  ~k~ r- 
ox 


:r=0 


(9) 


which  fixes  the  value  of  the  unknown  constant  A.  With  further  derivation,  omitted  for 
brevity,  it  can  be  shown  that  the  real  part  of  the  surface  temperature  is, 


m ,  .  In  LL  _  £  ,  X  7T 

r(l't)=  2 ke  "  cosM---4- 


(10) 


The  term, 


M  = 


(11) 


is  the  thermal  diffusion  length.  This  is  an  important  parameter  as  it  provides  a  length 
scale  for  the  penetration  depth  of  the  temperature  oscillations.  At  the  surface  of  the  solid, 


T(t)  =  cos(a ot - ). 

w  2k  K 


(12) 


2.1  Effect  of  Laser  Beam  Width 

The  preceding  theory  assumes  one-dimensional  (ID)  heat  flow,  which  is  generally  valid 
where  the  width  of  the  heat  source  is  large  compared  to  the  thermal  diffusion  length 
l_i.  This  condition  is  not  always  able  to  be  met  experimentally,  especially  when  laser 
irradiation  is  used  (the  most  common  type  of  source  used  for  TWI),  and  more  so  at 
low  modulation  frequencies  where  the  thermal  diffusion  length  can  be  relatively  large. 
Consider  for  example  the  inspection  of  an  aluminium  sample  at  a  modulation  frequency  of 
1  Hz.  Inserting  the  relevant  values  into  Eqn  11  shows  that  the  diffusion  length  in  this  case 
is  approximately  4  mm.  In  most  practical  situations  that  length  would  exceed  the  beam 
width,  resulting  in  strong  lateral  heat  flow  and  a  violation  of  a  key  assumption  of  Eqn  2. 

In  the  present  work,  the  influence  of  beam  diameter  on  the  phase  spectrum  was  exam¬ 
ined  by  means  of  numerical  simulation  using  the  finite  element  method  (FEM).  The  test 
case  involves  an  aluminium  slab  (a  =  52  x  10_5m2sec_1)  exposed  to  a  modulated  laser 
source  with  a  circularly-symmetric  Gaussian  beam  profile.  The  axisymmetric  geometry 
permits  a  2D  analysis.  At  each  modulation  frequency  the  simulation  time  was  adjusted 
to  allow  transients  in  the  response  to  extinguish  so  that  the  relative  phase  could  be  calcu¬ 
lated  from  the  steady-state  response.  Phase  spectra  were  computed  for  a  range  of  beam 
diameters.  Figure  3  compares  these  results  to  theoretical  predictions  derived  from  Eqn  2. 

The  influence  of  lateral  heat  flow  is  clearly  significant.  As  the  diffusion  length  in¬ 
creases  relative  to  the  beam  diameter  the  error  in  the  theoretical  phase  spectrum  grows, 
as  expected.  The  result  underscores  the  need  for  careful  experimental  design  if  analytical 
expressions  such  as  Eqn  2  are  to  be  used  as  a  basis  for  analysis.  These  numerical  findings 
are  consistent  with  experimental  results  reported  by  Fabbri  and  Cernuschi  [6].  Although 
lateral  heat  flow  is  an  impediment  to  analysis  it  does  not  diminish  TWI  as  a  material 
characterisation  tool.  Lateral  heat  flow  merely  introduces  an  additional  complexity  in  the 
phase  spectrum,  particularly  at  low  modulation  frequencies,  which  needs  to  be  factored 
into  analysis. 
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|i  (mm) 


Figure  3:  Phase  spectrum  as  a  function  of  laser  beam  diameter  computed  from  simulation 
(FEM).  Markers  trace  the  spectrum  corresponding  to  ID  heat  flow,  as  computed  from 
Eqn  2. 


3  Experimental  TWI  Facility 

Since  the  primary  objective  of  a  TWI  facility  is  an  accurate  measurement  of  the  relative 
phase  between  two  signals,  synchronisation  of  the  source  with  the  response  signal  is  critical. 
However,  any  real  measurement  inevitably  involves  time  delays,  which  are  a  potential 
source  of  bias  that  need  to  be  characterised  and  factored  into  the  computation  of  relative 
phase.  The  steps  taken  to  correct  for  such  factors  are  discussed  in  the  following. 


3.1  Apparatus  and  Equipment  Setup 

Figure  4  shows  a  photograph  of  the  experimental  setup.  The  main  components  of  the 
apparatus  are:  (i)  an  automated  micrometer-precision  x-y  stage  on  which  the  sample 
under  investigation  is  mounted,  (ii)  a  modulated  laser  source  for  irradiation  of  the  sample, 
and  (iii)  an  infrared  staring-array  camera  to  measure  the  induced  temperature  variations 
in  the  sample. 

Irradiation  of  the  sample  is  achieved  with  a  CUBE  660  nm  diode  laser  which  has  a 
maximum  output  power  of  100  rnW,  a  beam  divergence  of  approximately  1  milliradian, 
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Diode  Laser  Mirrors 


SC6000  Camera 


x-y  Stage 


Specimen 


Figure  4:  Experimental  apparatus  for  TWI. 


and  a  nominal  beam  diameter  of  4.8  mm1  measured  450  mm  from  the  source.  Preliminary 
tests  on  the  laser  revealed  that  the  beam  was  not  circularly  symmetric,  which  is  not 
unusual  for  diode  lasers.  Although  it  is  possible  to  reshape  the  beam,  the  process  involves 
a  trade-off  in  beam  power,  which  was  not  acceptable  for  the  present  work  where  the  laser 
had  to  be  operated  at  near  its  maximum  power  output  to  produce  a  satisfactory  response 
signal. 


The  infrared  camera  is  a  commercial  FLIR  SC6000  unit,  equipped  with  a  cryogenically 
cooled  640  x  512  staring-array  indium  antimonide  detector  with  a  sensitivity  of  approx¬ 
imately  17mK.  The  camera  produces  a  16  bit  digital  output  stream  which  is  acquired 
(nominally)  synchronously  with  a  reference  signal  using  custom  software  [7]  on  a  PC  run¬ 
ning  a  windows  operating  system.  The  reference  in  this  case  is  the  modulation  signal 
controlling  the  laser  output  which  was  supplied  by  a  stand-alone  function  generator. 


The  phase  difference  between  the  modulation  waveform  and  the  measured  temperature 
response  is  determined  from  the  argument  of  the  in-phase  and  quadrature  components 
of  the  cross-correlation.  The  correlation  sums  are  computed  in  real-time  which  results 
in  rapid  convergence.  Although  the  process  is  effective  in  rejecting  incoherent  signals 
(particularly  random  Gaussian  noise)  it  is  susceptible  to  systematic  phase  errors  caused 
by  time  delays  in  the  acquisition  system.  As  remarked  previously,  such  delays  have  the 
potential  to  bias  the  computed  phase  and  therefore  need  to  be  factored  into  the  cross¬ 
correlation  process. 


The  width  at  '2  of  the  maximum  power 
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3.1.1  Equipment  Characterisation 

The  output  data  stream  from  the  infrared  camera  includes  an  intrinsic  time  delay  asso¬ 
ciated  with  the  photon  integration  process,  as  well  as  a  delay  with  respect  to  the  skewed 
acquisition  of  the  reference  signal.  Essentially,  the  camera  signal  represents  a  photon  count 
over  a  prescribed  integration  time  r,  a  process  defined  by  the  following  frequency  response 
function, 


1 

T 


1 

LOT 


j'(cos(wr)  -  1) 


—  sin(cjr) 


The  phase  delay  is  then, 


(13) 


,  .  ,  /  1  —  cos(wr)\ 

4>t(lo)  =  arctan - ; — - — - —  .  (14) 

V  sin(uT)  ) 

In  the  SC6000,  the  integration  time  is  set  by  the  user,  which  means  that  (j)T(w)  is 
known.  Another  system  delay  results  from  a  latency  in  the  acquisition  of  the  infrared 
data  stream  and  the  reference  modulation  signal  which  are  acquired  by  separate  hardware 
components  in  the  PC.  This  latency  also  results  in  a  frequency  dependent  phase  offset, 
which  for  a  constant  delay,  can  be  expressed  as, 


4>td  =  ujS,  (15) 

where  td  is  an  acronym  for  “transport  delay”  where  that  delay  is  denoted  by  5  which  has 
the  unit  of  time.  Unlike  the  integration  time,  the  transport  delay  is  not  known  and  must 
be  determined  experimentally.  An  apparatus  described  previously  [7]  was  assembled  for 
that  purpose,  and  is  shown  schematically  in  Figure  5.  It  consists  of  an  area  black  body 
which  provides  a  uniform  and  stable  thermal  scene  that  is  transformed  to  a  harmonically 
modulated  scene  by  a  pair  of  linear  infrared  polarising  filters,  one  of  which  is  fixed  while  the 
other  rotates  at  a  prescribed  frequency  maintained  by  continuous  feedback  control.  The 
phase  of  the  radiation  transmitted  through  the  filter  assembly  is  determined  by  the  relative 
position  of  the  rotating  filter  which  is  tracked  by  a  photodiode/detector  unit  attached  to 
the  assembly,  as  shown  in  Figure  5.  The  arrangement  produces  an  arbitrary  phase  offset 
determined  by  the  angular  position  of  the  photodetector  unit  in  the  plane  of  rotation. 
Accordingly,  an  offset  term  4>0  is  added  to  the  phase  model,  which  becomes, 


4>  =  arctan( 


1  —  cos(a;T) 
sin(u;r) 


+  wS  +  (j)  0. 


(16) 


It  was  previously  mentioned  that  the  integration  time  is  set  by  the  user,  so  the  value 
is,  by  definition,  known.  However,  the  authors  did  not  verify  whether  the  prescribed 
value  was  precisely  implemented  by  the  camera.  Any  error  in  the  implementation  would 
invariably  flow  through  to  the  phase  model.  In  that  event,  r  would  need  to  be  treated  as 
unknown  and  evaluated  experimentally.  That  approach  was  not  taken  in  the  present  work 
where  the  prescribed  value  of  r  was  assumed  to  be  implemented  without  error.  With  that 
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assumption,  only  two  unknowns  need  to  be  determined  by  experiment:  <5  and  <f>0.  This 
was  achieved  by  means  of  linear  regression  applied  to  phase  measurements  adjusted  for 
the  integration  time,  which  amounts  to  minimising  the  following  objective  function, 


N  r* 
x  =  Y/Ui 

i=  1  V 

over  the  modulation  frequency 


—  arctan( - ^  1  ^ )  —  Wid  —  ,  (17) 

sm(o;jr)  J 

range  ui\  — >  lon  where  <fi  is  the  measured  phase  at  w*. 


Figure  5:  Experimental  setup  for  SC6000  camera  characterisation 


The  optimisation  was  done  over  a  discrete  set  of  frequencies  in  the  range:  1  —  100  Hz, 
in  increments  of  10  Hz.  The  measured  phase  at  each  frequency  was  derived  from  a  cross¬ 
correlation  applied  over  a  total  of  500  modulation  cycles.  Further  cross-correlation  pro¬ 
duced  no  additional  improvement  in  convergence.  Table  1  lists  5  and  (f)0  values  determined 
for  six  separate  experiments  conducted  under  precisely  the  same  conditions.  Both  values 
are  seen  to  vary  by  almost  10%.  The  large  variance  suggests  a  deficiency  in  the  phase 
model,  possibly  caused  by  the  assumption  of  a  constant  (time-invariant)  transport  delay, 
and  a  possible  problem  in  the  experimental  apparatus.  A  brief  investigation  of  the  vari¬ 
ability  was  made,  however  prior  to  discussing  its  outcome  it  is  worthwhile  illustrating  the 
effect  of  the  phase  corrections  defined  by  Eqn  16  as  implemented  for  the  average  values  of 
5  and  4>0-  Figure  6  plots  the  relative  phase  spectrum  measured  for  a  coated  metallic  sub¬ 
strate.  For  this  sample,  the  relative  phase  should  asymptote  to  —45°,  which  is  obviously 
not  achieved,  even  with  all  of  the  corrections  in  place,  however  the  corrections  are  seen  to 
significantly  reduce  the  error. 

A  quick  test  was  done  to  assess  whether  the  fluctuations  were  correlated  with  a  change 
in  ambient  temperature,  which  was  known  to  have  increased  monotonically  over  the  day. 
The  test  involved  measuring  the  transport  delay  and  offset  at  regular  intervals  across  a 
working  day.  The  values  are  shown  in  Table  2,  while  the  transport  delay  is  also  plotted  in 
Figure  7.  The  temperature  was  not  measured,  on  the  basis  that  the  variation  was  known 
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Figure  6:  Variation  of  phase  with  frequency  with  and  without  correction  for  system  delays. 
Table  1:  Transport  delay  and  phase  offset  measurements  from  initial  six  tests 


5  (ms) 

</>o  (°) 

-1.1047 

4.9143 

-1.0992 

4.7994 

-1.0923 

4.4641 

-1.0075 

4.6062 

-1.0099 

4.6034 

-1.0271 

4.8482 

to  be  relatively  large  («  5  °C)  and  monotonic  such  that  a  record  of  the  time  was  judged 
an  acceptable  surrogate.  The  results  in  Figure  7  provide  no  evidence  of  a  correlation. 

In  the  absence  of  a  model  for  the  error,  and  with  limited  time  available  for  further 
investigation,  a  pragmatic  solution  was  adopted  whereby  the  transport  delay  is  calculated 
as  part  of  the  phase  scanning  process.  The  approach  requires  a  reference  target  with  a 
known  phase  response  on  the  sample,  which  can  be  achieved  by  including  a  surface  layer 
with  a  thickness  that  is  some  factor  larger  than  the  thermal  diffusion  length,  such  that  the 
phase  difference  between  the  surface  response  and  the  source  modulation  approaches  45° 
(see  Equation  12).  Clearly,  a  check  is  required  to  ensure  the  beam  width  is  large  enough 
to  ensure  that  lateral  heat  flow  at  the  point  of  measurement  is  negligible. 

A  test  was  conducted  to  determine  the  extent  of  variation  in  the  system  transport 
delay  as  a  function  of  frequency.  A  metallic  sample  was  coated  in  a  layer  of  matt  black 
paint  ~  477 //m  thick2,  which  is  approximately  2.3  times  the  thermal  diffusion  length  at  a 
1  Hz  modulation  frequency  -  i.e.  the  coating  can  be  assumed  thermally  thick.  In  addition, 

2  as  measured  using  an  eddy  current  probe 
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Table  2:  Transport  delay  and  phase  offset  measurements  from  time  dependency  tests 


Time 

6  (ms) 

<fo  (°) 

09:22 

-1.0050 

4.4240 

09:59 

-1.0140 

4.7988 

10:42 

-0.9947 

4.5829 

11:13 

-0.9862 

4.5692 

11:49 

-1.0169 

4.0517 

13:59 

-1.0189 

4.0359 

14:35 

-0.9844 

4.2030 

15:10 

-1.0057 

4.8313 

15:46 

-1.0196 

4.0461 

16:30 

-1.0022 

4.5562 

£  -l 
<D 

Q 


^  -1.01 


_j _ i _ i _ i — 


09:22  09:59  10:42  11:13  11:49  13:59  14:35  15:10  15:46  16:30 

Time  (HH:MM) 


Figure  7;  Variation  in  SC6000  transport  delay  with  time  of  test. 


the  diffusion  length  is  small  relative  to  the  laser  beam  width  (~  4.8  mm)  so  the  assumption 
of  one  dimensional  heat  flow  should  be  well  met.  The  transport  delay  value  required  to 
achieve  a  relative  phase  measurement  of  —45°  was  recorded  for  modulation  frequencies  in 
the  range:  1  —  200  Hz.  These  values  are  shown  in  Figure  8.  The  result  raises  more  questions 
than  it  answers.  One  of  the  concerns  relates  to  the  strong  frequency  dependence.  This 
suggests  a  variable  latency  in  the  acquisition  of  the  separate  infrared  data  and  modulation 
waveform  data  streams  and  essentially  vindicates  the  in  situ  calibration  approach.  The 
second  concern  relates  to  the  fact  the  asymptotic  transport  delay  differs  significantly  from 
the  value  determined  using  the  calibration  apparatus.  The  reason  is  unclear  and  points 
to  the  need  for  further  work  on  addressing  this  issue.  Despite  these  unresolved  questions, 
it  was  decided  to  proceed  with  an  evaluation  of  the  system  in  measuring  film  thickness  on 
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Frequency  (Hz) 


Figure  8:  Variation  in  SC6000  transport  delay  with  frequency. 
a  prepared  test  sample. 


4  Paint  Thickness  Profiling 

Although  one  of  the  key  objectives  of  the  experimental  TWI  facility  is  to  furnish  a  basis  for 
the  measurement  of  adhesive  bondline  thickness,  that  particular  problem  involves  several 
factors  that  make  it  a  relatively  poor  test  case  for  a  preliminary  experimental  validation  of 
the  facility.  Instead,  a  slightly  simpler  test  case  was  conceived  that  has  important  similar¬ 
ities  to  the  original  problem.  The  case  entails  the  measurement  of  thickness  variations  in 
a  paint  coating  on  a  metallic  substrate.  It  is  attractive  for  the  following  reasons.  Firstly, 
as  a  surface  coating,  relatively  good  control  is  available  on  the  layer  thickness,  which  al¬ 
lows  it  to  be  closely  matched  to  the  bondline  thickness  in  a  typical  CMPLE  installation. 
Secondly,  independent  verification  of  the  layer  thickness  is  relatively  straightforward  since 
access  to  the  layer  is  unfettered.  Finally,  the  thermal  diffusivity  of  the  paint  chosen  for 
the  exercise  is  similar  to  that  of  epoxy,  ensuring  that  the  interference  behaviour  is  broadly 
representative  of  that  in  the  original  problem. 

4.1  Measurement  of  Paint  Thermal  Diffusivity 

A  measurement  of  paint  thickness  by  means  of  thermal  wave  interference  requires  an  a 
priori  knowledge  of  the  diffusivity.  Flash  Thermography  (FT)  was  used  to  determine 
the  thermal  diffusivity  of  the  paint  used  in  the  validation  exercise.  The  paint  is  an 
RS  matte  black  heat  radiator  paint  commonly  used  for  specimen  preparation  in  ther¬ 
mographic  nondestructive  inspection  and  thermoelastic  stress  analysis  applications.  A 
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through-transmission  configuration  was  adopted  for  consistency  with  Parker’s  [8]  flash 
method,  whereby  one  side  of  the  specimen  is  exposed  to  a  pulsed  heat  source  while  the 
rise  in  temperature  on  the  opposite  side  is  recorded.  In  order  to  produce  a  paint  film  that 
could  be  inspected  this  way,  a  carrier  was  prepared  consisting  of  a  rectangular  sheet  of 
cardboard  with  three  circular  apertures  each  covered  by  a  paint  film  of  different  thickness, 
as  shown  schematically  in  Figure  9.  The  films  were  created  by  spraying  paint  onto  a  thin 
layer  of  polyethylene  which  was  carefully  removed  from  the  carrier  once  the  paint  had 
dried.  The  thickness  of  each  layer  was  measured  using  an  eddy  current  probe.  Measure¬ 
ments  were  taken  at  several  positions  on  each  layer  and  averaged,  yielding  thickness  values 
of  45.85  ^rn(Li),  94.71  gm  (L2)  and  160.25  gm  (L3). 


L\  L2  L$ 

Figure  9:  Schematic  of  specimen  prepared  for  thermal  diffusivity  measurement. 


Parker’s  method  appears  straightforward,  however  it  assumes  a  thermal  impulse  re¬ 
sponse,  which  for  a  very  thin  object  is  difficult  in  practice  to  achieve.  The  problem  is 
illustrated  in  Figure  10,  which  compares  the  output  of  the  discharge  tube,  as  measured 
by  a  PbSe/Si  photo-detector,  with  the  rear-face  temperature  response  for  the  three  layers. 
The  duration  of  the  flash  output  is  not  negligibly  short,  although  a  reasonable  approxi¬ 
mation  is  achieved  for  the  thicker  paint  layer  L3. 

Parker’s  method  can  be  modified  to  accommodate  a  forced  response  by  convolving  the 
known  source  profile  with  the  impulse  response,  however  in  the  present  study  a  different 
approach  was  taken  that  produces,  in  effect,  the  same  result.  The  known  source  profile 
was  applied  in  an  FEM  simulation  to  generate  a  temperature  response  which  was  applied 
in  an  optimisation  process  where  the  thermal  diffusivity  was  varied  with  the  objective 
of  minimising  the  least-squares  error  between  the  measured  and  modelled  response.  The 
minimum  error  was  found  to  correspond  to  a  value  of  1.25  x  10“ '  m2s-1,  which  accords 
with  values  published  elsewhere  for  similar  paints  [9].  Figure  10  shows  excellent  agreement 
between  the  measured  response  and  simulations  based  on  that  value,  particularly  in  the 
early  stages  of  the  response  when  through-thickness  diffusion  is  the  dominant  heat-transfer 
mechanism.  Note  that  the  peak  in  the  temperature  response  marks  a  switch  from  conduc¬ 
tion  to  convection  and  radiation  dominated  heat  transfer,  which  explains  the  discrepancy 
that  develops  between  the  measured  and  modelled  response  for  the  thin  sample  as  time 
progresses.  The  objective  function  was  windowed  to  exclude  that  portion  of  the  response 
so  this  discrepancy  does  not  effect  the  estimate  obtained  for  a. 
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Figure  10:  Flash  profile,  temperature  traces  recorded  for  L\,  L2  and  L3;  as  well  as  cor¬ 
responding  FEM  simulations  for  the  optimal  value  of  a. 


An  estimate  for  the  thermal  diffusivity  was  also  derived  using  Parker’s  method  by 

inserting  the  half-maximum  rise-time  (denoted  by  1 1)  for  the  thin  layer  response  into  the 

2 

expression  [8], 


a  = 


1.38L2 

TT2t  1 


(18) 


An  estimate  of  8.9  x  10-8  m2s-1  is  obtained,  which  differs  from  the  optimised  value 
by  almost  30%,  reinforcing  the  need  to  consider  the  duration  of  the  flash  in  the  diffusivity 
calculation. 


4.2  Non-contact  Measurement  of  Paint  Thickness 

To  facilitate  an  experimental  evaluation  of  the  performance  of  the  facility  in  characterising 
film  thickness,  a  test  specimen  was  produced  with  a  coating  of  the  matte  black  radiator 
paint  described  in  the  previous  section.  The  substrate  was  a  1.6  mm  thick  rectangular 
coupon  of  A12024,  roughly  150  x  40nun  in  size.  A  stepped  coating  thickness  profile  was 
produced  by  applying  the  paint  as  an  aerosol  in  multiple  passes  where  strips  approximately 
30  mm  wide  were  masked  once  the  desired  number  of  paint  layers  had  been  applied.  The 
end  result  is  shown  schematically  in  Figure  11.  An  eddy  current  probe  was  used  to  measure 
the  coating  thickness  for  each  strip.  Probe  measurements  were  taken  at  several  randomly 
selected  points  within  each  strip  and  a  mean  and  variance  calculated,  which  are  listed 
in  Table  3.  Note  that  the  sixth  strip  was  made  substantially  thicker  in  order  to  furnish 
the  “half-plane”  reference  target  required  for  the  transport  delay  calibration  described 
previously. 
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Li 


L2  L3  L4 


L5 


Lq 


Figure  11:  Plan  and  front  views  of  Al202f  aluminium  alloy  specimen  with  six  different 
paint  coating  thicknesses  (dimensions  not  to  scale). 


Table  3:  Thickness  measurements  in  microns  taken  with  an  eddy  current  probe. 


Thickness 

Li 

L2 

L3 

l4 

L5 

La 

h 

20.44 

44.12 

74.87 

88.07 

100.52 

476.84 

a 

0.97 

1.63 

2.07 

2.64 

3.33 

7.42 

p  —  a 

19.47 

42.49 

72.80 

85.43 

97.19 

469.42 

g  +  cr 

21.41 

45.75 

76.94 

2.64 

90.71 

484.26 

Interference  occurs  only  when  the  diffusion  length  is  of  the  same  order  as  the  film 
thickness  which  means  that  the  performance  of  TWI  varies  significantly  with  modulation- 
frequency.  Since  the  material  properties  of  the  layered  sample  are  known,  as  is  the  thick¬ 
ness  profile,  at  least  approximately,  Eqn  2  can  be  used  to  determine  a  priori  a  set  of 
optimal  frequencies  for  inspection.  These  frequencies  were  determined  by  calculating  the 
phase  spectrum  for  a  layer  thickness  one  standard  deviation  above  and  one  standard  de¬ 
viation  below  the  mean  measured  paint  thickness.  The  two  spectra  are  expected  to  differ 
most  at  precisely  the  optimal  inspection  frequency.  That  approach  was  applied  to  strips 
Li,  L3  and  L5,  yielding  frequencies  of  4,  7  and  91.2  Hz  respectively  (Figure  12). 

The  photograph  in  Figure  4  shows  the  sample  fixed  to  the  x-y  stage  in  preparation  for 
a  scan.  A  phase  profile  was  formed  by  measuring  the  phase  at  discrete  points  spaced  2  mm 
apart  along  the  scan  line  shown  in  Figure  11.  Each  discrete  measurement  was  formed  from 
a  cross-correlation  involving  between  250  and  1000  modulation  cycles,  depending  on  the 
modulation  frequency  -  in  general,  higher  frequencies  requiring  more  cycles.  Scans  were 
performed  at  the  three  optimal  frequencies,  as  well  as  at  an  arbitrary  frequency  of  30  Hz, 
which  was  added  for  comparative  purposes. 

The  measured  phase  profiles  are  shown  in  Figure  13.  Overall,  the  phase  variation  across 
the  sample  is  significant,  with  a  strong  discontinuity  in  phase  between  each  strip.  The 
profile  suggests  an  uneven  thickness  within  some  of  the  strips,  which  is  most  evident  for  L3. 
The  size  of  the  variation  for  L3  was  not  anticipated  given  the  care  taken  in  preparing  the 
specimen.  It  is,  however,  encouraging  to  note  that  these  variations  are  largely  consistent 
in  form  across  the  four  frequencies  but  amplified  for  the  frequencies  predicted  by  theory 
to  provide  optimal  sensitivity  at  the  nominal  layer  thickness  -  e.g.  the  variation  in  phase 
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Figure  12:  Calculated  phase  spectra  (solid)  and  their  d\ 
timal  modulation  frequencies  for  three  different  paint  thii 
and  Lk  (bottom). 
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across  strip  L3  is  greatest  for  a  modulation  frequency  of  7  Hz. 

Phase  profiles  corresponding  to  the  thickness  measurements  listed  in  Table  3  were 
calculated  using  FEM  and  Eqn  2.  These  synthetic  profiles  are  plotted  alongside  the 
measured  profiles  in  Figure  14.  The  agreement  between  the  three  sets  of  results  is  generally 
good,  noting  that  the  synthetic  profiles  correspond  to  an  average  strip  thickness  so  cannot 
describe  the  thickness  variations  within  each  strip.  Notwithstanding  this  limitation  some 
important  subtleties  are  captured.  Note  for  instance  the  91.2  Hz  profile  and  the  small 
increase  in  phase  from  L2  to  L3.  The  increase  is  considered  significant  as  it  runs  against 
the  overall  trend  of  a  monotonically  decreasing  phase  at  the  lower  modulation  frequencies. 
A  similar  small  phase  increase  against  trend  is  seen  for  the  transition  between  L4  and  L5 
for  the  30  Hz  result. 


Figure  13:  Phase  measurements  across  the  coated  sample. 


A  more  rigorous  validation  of  the  experimental  result  requires  that  the  phase  profile 
produced  by  the  facility  be  converted  to  a  thickness  profile,  since  this  then  allows  a  com¬ 
parison  to  be  made  against  independent  measurements  of  thickness.  The  two  profiles  are 
related  through  heat  diffusion  so  the  conversion  requires,  in  effect,  the  solution  of  an  in¬ 
verse  mathematical  problem,  where  the  model  can  be  furnished  by  Eqn  2  if  the  heat  flow 
is  largely  one-dimensional,  or  by  FEM  where  lateral  heat  flow  is  significant.  The  latter 
approach  was  used  in  the  present  study.  A  Nelder-Mead  Simplex  optimisation  routine 
applied  in  Matlab  was  used  to  solve  the  relevant  inverse  problem.  Four  thickness  profiles 
were  produced,  corresponding  to  the  four  sets  of  experimental  data.  The  profile  estimates 
are  shown  in  Figure  15.  The  thickness  profile  is  of  course  unique  so  the  ambiguity  in 
these  estimates  is  evidence  of  error,  relating  either  to  the  inversion  process  and  more  likely 
deficiencies  in  the  experimental  data.  Some  of  the  error  is  easily  explained.  For  example, 
the  absence  of  a  profile  in  the  91.2  Hz  result  probably  reflects  weak  interference  caused 
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91.2  Hz  30  Hz 


7  Hz  4  Hz 


Figure  14:  Comparison  between  measured  (markers),  FEM  (dashed  line)  and  theoretical 
(solid  line)  phase  profiles. 


by  the  short  diffusion  length  at  this  frequency  relative  to  the  layer  thickness,  which  is 
supported  by  the  largely  consistent  behaviour  in  the  measured  phase  profile  within  the 
same  region.  The  discrepancy  in  the  30  Hz  result  is  similarly  explained,  noting  that  the 
erroneous  behaviour  occurs  over  a  smaller  range  of  coating  thickness  relative  to  the  91.2  Hz 
result,  consistent  with  deeper  penetration  at  the  lower  frequency.  Elsewhere,  the  discrep¬ 
ancy  between  profiles  is  much  smaller  and  mostly  characterised  by  a  relatively  small  offset 
in  the  profile  which  is  otherwise  consistent. 

As  none  of  the  modulation  frequencies  furnish  a  valid  thickness  profile  across  the  entire 
sample  it  was  decided  to  assemble  a  complete  profile  from  the  results  corresponding  to  the 
frequency  optimised  for  the  particular  region,  i.e.  the  profile  at  91.2  Hz  was  selected  for 
Li,  the  profile  at  7 Hz  was  used  for  L3  and  the  profile  at  4 Hz  was  chosen  for  L5.  As  the 
nominal  thicknesses  of  L4  and  L5  are  similar,  the  profile  at  4  Hz  was  also  used  for  L4.  The 
profile  at  30  Hz  was  used  for  L2  on  the  basis  that  that  frequency  sits  between  the  optimal 
frequencies  determined  for  the  strips  to  either  side. 

A  further  independent  measurement  of  the  paint  thickness  profile  was  obtained  using 


UNCLASSIFIED 


17 


DSTO-TR  2760 


UNCLASSIFIED 


Figure  15:  Paint  thickness  profiles  generated  by  FEM  for  each  of  the  four  frequencies. 


a  surface  texture  measuring  instrument3,  which  we  refer  to  now  on  as  a  profilometer.  The 
system  uses  a  tracing  arm  attached  to  the  drive  unit  housing  via  a  magnetic  bearing.  Ac¬ 
cording  to  the  technical  specifications  of  the  instrument,  it  yields  a  surface  measurement 
accuracy  of  0.04  /jrn.  Profiles  were  recorded  along  three  parallel  scan  lines,  one  corre¬ 
sponding  approximately  to  the  phase-profile  scan  line  shown  in  Figure  11,  and  another 
two  taken  along  lines  1  mm  above  and  1  mm  below  the  original.  These  additional  scans 
were  done  to  compensate  for  uncertainties  in  the  position  of  the  phase-profile  scan  line. 

Figure  16  compares  these  independently  sourced  profiles  to  the  TWI  result  and  the 
average  strip  thickness  measured  using  the  eddy  current  probe.  Part  of  the  profile  for  L5  is 
omitted  as  the  surface  texture  measuring  instrument  was  limited  in  its  traversal  to  120nun. 
The  agreement  overall  is  excellent.  An  important  observation  here  is  that  the  intra-strip 
variations  in  the  TWI  profile  also  occur  in  the  profilometer  scans,  proving  that  these  are 
variations  in  thickness  rather  than  artefacts  of  the  measurement  system.  In  making  that 
assessment  allowance  needs  to  be  made  for  the  large  scatter  in  the  profilometer  scans, 
which  is  thought  to  be  noise  rather  than  a  physical  feature  of  the  coating.  Analysis  of  the 
scatter  shows  that  it  accounts  for  a  standard  deviation  of  approximately  1.5  gm  which  is 
an  order  of  magnitude  above  the  nominal  accuracy  of  the  device.  By  comparison,  the  TWI 
profile  is  much  cleaner,  suggesting  a  noise  floor  significantly  better  than  1.5  //rn.  There  are 
also  some  concerning  aspects  about  the  performance  of  TWI.  Within  the  thinnest  strip, 
L] ,  TWI  furnishes  a  thickness  estimate  20%  lower  than  the  profilometer  and  eddy  current 
probe,  which  is  a  significant  discrepancy  that  implicates  the  system  phase  error  discussed 
previously. 


3Mahr  MarSurf  XC10  system  with  a  MarSurf  CD  120  drive  unit 
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Figure  1 6:  Comparison  of  thickness  profiles  derived  from  an  eddy  current  probe,  a  TWI 
profile  estimate,  and  a  surface  texture  measuring  instrument  (three  separate  scans). 

5  General  Discussion 

An  error  in  the  transport  delay  is  known  to  manifest  an  error  in  the  system  phase  output 
that  grows  with  frequency  (see  Eqn  15),  which  was  illustrated  in  Figure  6  where  the  phase 
error  was  obvious  for  frequencies  above  200  Hz.  Unfortunately,  it  was  not  possible  to 
conclusively  determine  the  phase  error  within  the  range  of  frequencies  considered  in  the 
test  case,  however  the  fact  the  greatest  discrepancy  in  Figure  16  occurs  for  the  highest 
modulation  frequency,  i.e.  91.2  Hz,  is  consistent  with  a  transport  delay  error. 

There  is  reason  to  also  question  the  fidelity  of  the  FEM  model  used  in  translating  the 
measured  phase  values  to  a  coating  thickness.  In  the  model,  the  beam  was  assumed  to  have 
a  circularly  symmetric  Gaussian  intensity  profile  which  is  known  to  be  an  idealisation  of 
the  true  profile.  Figure  17  provides  an  approximate  measurement  of  the  true  profile  formed 
by  mapping  the  relative  phase  across  the  beam  on  a  thermally  thick  sample4.  A  contour 
is  included  to  clarify  the  asymmetry  of  the  profile.  The  beam  distortion  means  that  the 
model  is  likely  to  have  understated  the  impact  of  lateral  heat  flow.  The  associated  phase 
error  should  increase  as  the  modulation  frequency  decreases,  as  illustrated  in  Figure  3, 
so  if  the  effect  has  had  a  significant  bearing  on  the  results  in  Figure  16  it  should  be 
most  apparent  for  the  thicker  coatings.  Such  behaviour  was  not  observed,  suggesting  the 
discrepancy  at  L\  is  unlikely  to  relate  to  lateral  heat  flow  and  is  more  probably  associated 
with  uncertainty  in  the  system  transport  delay,  as  described  previously. 

In  view  of  the  obvious  shortfalls  in  system  performance,  the  results  in  Figure  16  provide 
a  conservative  guide  to  the  full  potential  of  the  technique.  Given  the  remarkable  quality 
of  the  result  the  potential  is  evidently  large,  and  gives  considerable  cause  for  optimism  in 
achieving  the  original  aim  of  the  facility,  which  is  to  provide  an  experimental  platform  for 
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Figure  17:  Asymmetry  in  the  phase  across  irradiated  spot  on  a  thermally  thick  sample. 


the  development  of  an  in  situ  bondline  thickness  sensor. 


6  Conclusions 


This  report  has  detailed  the  development  of  a  thermal  wave  interferometry  system  for  the 
characterisation  of  film  thickness.  The  system  furnishes  a  measurement  of  the  relative 
phase  between  a  modulated  heat  flux  applied  to  a  thin  film  by  a  diode  laser  and  the 
oscillatory  temperature  response  of  the  film.  The  phase  measurement  is  converted  to  a 
thickness  estimate  by  means  of  an  inverse  solution  of  the  heat  diffusion  equation.  An 
accurate  measurement  of  the  phase  proved  difficult  to  achieve  as  latencies  in  the  data 
acquisition  system  were  found  to  vary  with  modulation  frequency,  as  well  as  randomly. 
Corrections  for  these  latencies  were  developed  and  implemented  and  proved  successful  in 
part. 

A  validation  of  the  facility  was  performed  on  an  experimental  test  coupon  prepared 
with  paint  layers  of  varying  thickness.  Estimates  of  the  paint  thickness  profile  furnished 
by  the  facility  were  compared  to  profiles  measured  by  an  eddy  current  probe  and  a  surface 
profilometer.  Both  provided  corroboration  for  key  aspects  of  the  TWI  result,  however 
evidence  was  also  produced  suggesting  deficiencies  in  the  approach  used  to  compensate 
for  transport  delays  in  the  acquisition  system.  In  view  of  the  considerable  scope  for  im¬ 
provement  in  the  phase  stability  of  the  system,  the  fidelity  of  the  result  was  considered 
remarkable.  Work  is  currently  underway  on  making  these  improvements,  with  the  antic¬ 
ipation  that  they  should  result  in  a  film  thickness  measurement  threshold  of  well  below 
1  ^m.  Once  that  is  achieved  work  will  proceed  on  the  development  of  an  in  situ  TWI 
sensor  that  can  be  incorporated  into  bonded  strain-sensing  devices  such  as  the  CMPLE. 
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